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A coupled-channel model previously employed to describe the narrow _D*o(2317) and broad 
-Do (2400) charmed scalar mesons is generalized so as to include all ground-state pseudoscalar- 
pseudoscalar and vector-vector two-meson channels. All parameters are chosen fixed at published 
values, except for the overall coupling constant, which is fine-tuned to reproduce the _D*o(2317) 
mass. Thus, the radial excitations _Dso(2850) and _Do(2740) are predicted, both with a width of 
about 50 MeV. The former state appears to correspond to the new _D3,j(2860) resonance decaying 
to DK announced by BABAR in the course of this work. Also the -Do(2400) resonance is roughly 
reproduced, though perhaps with a somewhat too low central resonance peak. 

PACS numbers: 14.40.Lb, 14.40.Ev, 13.25.-k, 12.39.Pn 



The discovery of the L>;o(2317) 0] and Dsi{2A60) ^ 
charm-strange mesons three years ago has triggered a 
strongly renewed interest in heavy-light mesons, and even 
meson spectroscopy in general. Especially the £'*q(2317) 
has given rise to many different theoretical efforts (see 
Ref. I3] for a long though still incomplete list of refer- 
ences). The reason is its surprisingly low mass, some 
170 MeV below the predictions of standard relativized 
constituent quark models for the ground-state scalar cs 
meson (see e.g. Ref. 0|), an assignment that has in the 
meantime been confirmed by experiment Q. This dis- 
crepancy led several model builders to propose alterna- 
tive explanations for the £'*g(2317), such as a tetraquark 
or a meson molecule. However, in Ref. 01 we showed 
how the low mass of the £'*o(2317) can be quantitatively 
understood by taking into account its strong coupling to 
the nearby S-wave DK channel. This explanation was 
later supported by Refs. and 0. Similarly, we ex- 
plained the Dsi (2460) in Ref. via its strong coupling to 
the S'-wave D*K threshold. The coupled-channel model 
employed in Ref. had been previously used, with es- 
sentially the same parameters, to reproduce the S'-wave 
Kn phase shifts and predict [l^ the K^iSOO) (alias k) 
resonance, later confirmed by experiment [5j- 

Nevertheless, no consensus has been reached so far on 
the D*o(2317), in part due to the poor experimental sta- 
tus of the very broad partner charm-nonstrange state, 
listed as -Do(2400), but first reported at a mass of 2308 
MeV ^nj and later also at 2407 MeV Therefore, a 
more detailed coupled-channel analysis of charmed scalar 
mesons is very opportune, also in view of new and heav- 
ier states that are expected to be found at B factories. 
Clearly, for a reliable description of higher resonances, 
additional decay channels must be accounted for. Thus, 
in the present Letter, we extend 0| the model of Ref. ^ 



by including all lowest pseudoscalar-pseudoscalar (PP) 
and vector-vector (VV) two-meson channels that couple 
to the scalar cs and cn (n = u, d) systems, in an ap- 
proach very similar to Ref. In the latter paper, the 
coupling to all PP channels allowed to fit the properties 
of the light scalar mesons cr, k, ao(980), and /o(980), such 
as phase shifts, line-shapes, elasticities, and inelastic am- 
plitudes, obtaining an overall good description of these 
observables, as well as very reasonable pole positions. In 
the present investigation, the inclusion of the VV chan- 
nels as well is crucial to study possible radial excitations, 
as in the cn and cs sectors the lowest VV channels open 
at roughly 2.8 GeV and 2.9 GeV, respectively. 

We will first discuss the results and finish this Let- 
ter with a short description of the mathematics behind 
the Resonance-Spectrum Expansion (RSE), which is the 
framework of our model 15]. In Fig. ^ we show, for the 
cs case, the resulting S-wave DK DK cross sections. 
The dashed line refers to the case where only PP chan- 
nels are included, the solid line to the case where also 
VV channels are accounted for. We will discuss the latter 
case. A comparison of differences for the two situations 
is presented below for the cn system. 

At energies close to threshold (at 2.363 GeV), the cross 
sections are large due to the presence of the 15*0(2317) 
bound state just below threshold. For higher total invari- 
ant mass {\/s), the cross sections decrease, however not 
as fast as expected, due to the presence of a scattering- 
matrix pole, which we find at 2779- i233 MeV. At 2.516 
GeV one observes the effect of the opening of the Dsf] 
channel, while at about 2.79 GeV the cross sections al- 
most vanish. The first radial excitation of the D*o(2317) 
is found with a peak mass of 2847 MeV and a width of 
47 MeV, and so i s a g ood candidate for the new BABAR 
state Z?sj(2860) '16'|, which decays to DK and not to 
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FIG. 1; The predicted S-wave DK DK cross sections. 
The dashed curve corresponds to PP channels only, the solid 
curve to PP + VV channels. In the inset, we compare the 
model results to the data, with arbitrary normalization taken 
from Ref. Q. 

D*K, having a mass of 2857 MeV and a width of 48 
MeV. In our model it is associated with a resonance pole 
at 2842 - 123.6 MeV. From the inset of Fig. [2 one can 
judge how well our D*o'(2850) predicts the line shape of 
BABAR's Ds j{2860). There is furthermore some indica- 
tion .16j that the data need a broad state as well, which 
might correspond to our pole at 2779 — i233. 

In Fig. 121 we show, for the cn case, the resulting S- 
wave Dtt Dtt cross sections. We find the lowest reso- 
nance pole at 2149-ailMeV (PP) or 2174-i96.4 MeV 
(PP-hVV), with peak mass at 2180 MeV or 2190 MeV, 
respectively. This broad resonance should correspond 
to the £)g(2400). Our prediction seems too low, but is 
not unreasonable in view of the unsettled experimental 
situation Q, and also considering our highly dynami- 
cal Z?q(2400) pole ,6j, which can travel a long distance 
with moderate changes in the model's coupling constant 
(A). For instance, if we reduce A somewhat so as to let 
the Z)*g(2317) become slightly heavier, though still below 
the DK threshold, it is possible to increase the Dq(2400) 
mass prediction by up to 100 MeV. Nevertheless, we 
beheve it is safer to keep the established I?*q(2317) in 
its place, considering the persisting uncertainties regard- 
ing the Z?g(2400). Moreover, the experimental values 
concern production processes, and not elastic scattering. 
Furthermore, the analyses rely on a Breit-Wigner shape 
for the 1)2(2460) resonance, which has a large contribu- 



tion to the total signal. In the inset of Fig. El we show a 
comparison of our signal with the data |ll| , for invariant 
masses well below the -D2(2460) resonance. 

Next we look for poles at higher energies. In this sit- 
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FIG. 2: The predicted elastic S-wave Dtt cross sections. The 
dashed curve correspond to PP channels only, the solid curve 
to PP -I- VV channels. In the inset, we compare the model 
results to the data, with arbitrary normalization taken from 
Ref. m. 

nation, we only consider the poles for the full PP+VV 
system, as several VV channels open above ~ 2.8 GeV. 
Still in the cn case, we find a relatively narrow pole at 
2737 - ■j24.0 MeV and a very broad one at 2703 - i228 
MeV. The narrow state, with a width of about 50 MeV, 
corresponds to the first radial excitation of the cn system, 
shifted to complex energy by the coupled channels, while 
the very broad resonance is the strongly distorted and 
shifted ground state of the confinement spectrum, also 
found in Ref. , though now with a width of roughly 450 
MeV instead of « 200 MeV. Note that the (2400) is a 
dynamical continuum pole, just as in Ref. The nar- 
row resonance at about 2.74 GeV predicted here should 
be observable, though the S'-wave elastic Dtt cross section 
is quite small (see Fig. (21 solid curve) . 

Finally, let us turn to a short description of the model 
employed in this work. The model's scattering matrix 
(5) for N two-meson channels (masses Mn and Mi2, i ~ 
1, 2, 3, . . ., N, orbital angular momentum £i, relative 
linear momentum ki), all coupled to one quark- antiquark 
confinement channel (radial confinement spectrum given 
hy Eo = uj {£gg + 3/2) + niq + mq, Ei = Eq + 2u, E2 = 
El + 2u!, . . .), has the following closed form [l5l |: 
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1 + 2iX^rqg ^ ^ < ^ ^ — — — ? fJ'mkmjim {kmTqq) {kmTqq) 



where, in the i-th channel, the hnear momentum ki and 
reduced mass /Xi are related to the total invariant mass 
E of the system, and to the two meson masses Mn and 
Mi2, through 

E = ^Jkf + Mf, + ^kf + Mf, , 
= 2kf + Ml+M^^ + 2Efi, . 

ji and /i^^'' represent the spherical Bessel and Hankel 
function of the first kind, respectively. 

The model parameters representing quark masses 
(m„ = 0.406 GeV, = 0.508 GeV, ruc = 1.562 GeV) 
and the radial spacings in the bare confinement spectrum 
{uj = 0.19 GeV) are kept identical to the ones originally 
optimized in Ref. ^^^o used in Ref. More- 

over, the parameter r^g, which stands for the average ra- 
dius of ^Pq quark-pair creation, is identical to the value 
= 3.2 GeV^^ used in Ref. 10], but scaled with the 
reduced quark mass in order to impose flavor symmetry 
of our equations 

TOc (TOs + rrin ) 

= Uhlpllll^r^.^ . 1.88 GeV- . 
nic [mn + ms) 

The overall decay coupling constant A is fine-tuned to 
reproduce the mass of the now very well established 
Z?*g(2317). Yet, also A turns out to be close to the values 
used in the light scalar sector 0| , owing to the referred 
flavor-symmetric mass scaling. This yields the values 
A = 2.854 GeV-3/2 ^^len only PP channels are included, 
and A = 2.617 GeV-3/2 with PP as weU as VV channels. 
Note that the former value of A is fully compatible with 
the values found for the light scalars in Ref. fl3|, which 
analysis was also restricted to PP channels. The change 
in A from the VV channels amounts to a reduction by 
less than 10%. 

The channels included in the present work are summa- 
rized in Tabled their relative couplings to the qq channels 
in Table in 

The pseudoscalar 77-77' mixing angle Ops we choose at 
the recently found experimental value Qps — —13.5° jl^l 
(octet-singlet basis). However, we also verify our results 
for another frequently used value, i.e., Qps = —17.3° 
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Channels 


Thresh. 


(waves) 


(GeV) 


(waves) 


(GeV) 


Dn (S) 


2.004 


DK (S) 


2.363 


Dv (S) 


2.415 


Dsn {S) 


2.516 


Dv' (S) 


2.825 


Dsv' (S) 


2.926 


DsK {S) 


2.464 






D-p (S,D) 


2.784 


D*K* {S,D) 


2.902 


D'uj {S,D) 


2.791 


D*s<j} {S,D) 


3.132 


DIK' {S,D) 


3.006 







TABLE I: The various meson-meson channels included in this 
analysis, and their threshold energies. 



charm-nonstrange 


charm-strange 


v/1/16 


v/1/12 


x^l/144-y^l/72 


Vl/72 - a;^l/36 


^^1/144 + ^^1/72 


x^l/72-yVl/36 


v/l/24 




v/l/48 , v/5/12 


v/l/36 , v/579 


Vl/144 , v/5/36 


v/l/72,^5/18 


Vl/72 , v/5/18 





TABLE II: The relative couplings (four-fermion recombina- 
tion coefficients 1'^'' various meson-meson channels 
included in this analysis (in the same order as in Table m, 
to J^*^ = 0^"'' cn and cs in 5* or _D waves. For the rela- 
tive couplings to higher radial excitations n, one has g{n) = 
y^{n + l)g{0)/2" for S and g{n) = y/{2n + 5)g{0)/{5 x 2") 
for D waves. The symbols x and y stand for cos 6pg and 
sinOpg, respectively. 

|l4l |. which turns out to change the predictions by only a 
few MeV. We force the damping of closed scattering chan- 
nels with subthreshold form factors, which are a standard 
tool in modern multichannel phase-shift analyses: 

^ 5?(")e"^'^ for ^ekf<0 . (2) 
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We choose the value a = 4 GeV~^, which is the same 
as used in the analysis of the light scalars Such 
a suppression, in addition to the one resulting from our 
kincmatically relativistic Schrodinger formalism, can be 
justified from relativistic covariance, offshellness, self- 
energies, and other effects not accounted for in the 
present model. These contributions are, of course, very 
difficult to rigorously evaluate in our nonperturbative 
scheme. However, even if we were to completely switch 
off subthreshold damping, our Z?*q(2850) pole would only 
shift to 2864 — i x 15 MeV, after a readjustment of A so 
as to reproduce again the £)*g(2317) mass. 

Having now fixed all parameters in formula Q , we can 
search our amplitudes for resonance poles. We predict 
the first radial excitations of the D*q{2317) and D^{2A00) 
to come out as D*J{2850) ^ and I?S'(2740), respec- 
tively. We furthermore predict the very broad states 
i:»;o(2780) and D*{2700), which might show up in a more 
pronounced way in production experiments than in elas- 
tic scattering. 

Note added in proof: after completion of this work, 
the BABAR collaboration posted a preprint HJ] con- 
firming the announcement of the Dsj(2860) in Ref. 
Furthermore, three theoretical papers on this new state 
have appeared in the meantime, the first one |22l | favor- 
ing a 3- (1 ^Da) assignmentjhe second (H a 0+ i2^Po) 
like we do, and the third j2^ admitting either possibil- 
ity. Clearly, the non-observation so far of the D*K decay 
mode, forbidden for a scalar meson, favors the 0^ option, 
although the DsT] mode, not observed either, is allowed in 
both the 0+ and 3~ scenarios. It is also interesting that 
Ref. jHI, which makes out a case for the 3~ assignment, 
predicts branching ratios Ds ji2860) D*K/DK = 
0.39 for 3-, and Dsj(2860) DsTj/DK = 0.34 for 0+. 
We predict a value of 0.30 for the latter branching ratio, 
if we include all PP+VV channels. Anyhow, experiment 
will have the final word on interpreting the Dgj (2860) 
beyond any doubt, by observing either Dst] or D*K. 



Acknowledgments 

This work was presented at the QNP06 conference in 
Madrid |l3l | , unknowing of the simultaneous experimental 
release at the CHARM06 conference in Beijing We 
are indebted to S. Tosi for drawing our attention to the 
brand-new data. This work was supported in part by the 
Fundagao para a Ciencia e a Tecnologia of the Ministerio 
da Ciencia, Tecnologia e Ensino Superior of Portugal, 
under contract POCI/FP/63437/2005. 



[1 
[2: 

[s: 
[4: 

[6; 

[7: 

[8 

[9: 
[10: 

[11 

[12 
[13 



[14 

[is: 
[16: 

[17 

[is: 

[19 
[20 

[21 
[22: 
[23: 

[24: 

[25: 



B. Aubert et al. [BABAR Collaboration], Phys. Rev. 
Lett. 90, 242001 (2003) arXiv:hep-ex/0304021 . 

D. Besson et al. [CLEO CollaborationJ , Phys. Rev. D 68, 
032002 (2003) arXiv:hep-ex/0305100 . 

E. van Beveren, J. E. G. N. Costa, F. Kleefeld, 
and G. Rupp, Phys. Rev. D 74, 037501 (2006) 
arXiv:hep-ph/0509351 . 

S. Godfrey and N. Isgur, Phys. Rev. D 32, 189 (1985). 
W. M. Yao et al. [Particle Data Group], J. Phys. G 33, 
1 (2006). 

E. van Beveren and G. Rupp, Phys. Rev. Lett. 91, 012003 

(2003) arXiv:hep-ph/0305035 . 

D. S. Hwang and D. W. Kim, Phys. Lett. B 601, 137 

(2004) arXiv:hep-ph/0408154 . 

Yu. A. Simonov and J. A. Tjon, Phys. Rev. D 70, 114013 
(2004) arXiv:hep-ph/0409361 . 

E. van Beveren and G. Rupp, Eur. Phys. J. C 32, 493 

(2004) arXiv:hep-ph/0306051 . 

Eef van Beveren and George Rupp, Eur. Phys. J. C 22, 
493 (2001) arXiv:hep-ex/0106077 . 

K. Abe et al. [Belle Collaboration], Phys. Rev. D 69, 

112002 (2004) arXiv:hep-ex/0307021 . 

J. M. Link et al. [FOCUS Collaboration], Phys. Lett. B 

586, 11 (2004) arXiv:hep-ex/0312060 . 

E. van Beveren and G. Rupp, Multichannel calculation of 

the very narrow D* j (2317) and the very broad Dq (2300- 

2400), parallel talk given by G. Rupp on 6 June 2006 at 

the IVth Intern. Conf. on Quarks and Nuclear Physics, 

5-10 June 2006, Madrid, Spain. 

E. van Beveren, D. V. Bugg, F. Kleefeld and G. Rupp, 
Phys. Lett. B 641, 265 (2006) arXiv:hep-ph/0606022 . 
E. van Beveren and G. Rupp, Int. J. Theor. 
Phys. Group Theor. Nonlin. Opt. 11, 179 (2006) 
arXiv;hep-ph/0304105 . 

A. Palano, New Spectroscopy with Charm quarks at B 
factories, plenary talk given on 7 June 2006 at the 
Charm2006 Intern. Workshop, 5-7 June 2006, Beijing, 
China; Silvano Tosi, private communication. 
E. van Beveren, G. Rupp, T. A. Rijken, and C. DuUe- 
mond, Phys. Rev. D 27, 1527 (1983). 

E. van Beveren and G. Rupp, Mod. Phys. Lett. A 
19, 1949 (2004) arXiv:hep-ph/0406242 ; G. Rupp, 

F. Kleefeld, and E. van Beveren, AIP Conf. Proc. 756 

(2005) 360 arXiv:hep-ph/0412078 ; F. Kleefeld, AIP 
Conf. Proc. 717, 332 (2004) arXiv:hep-ph/0310320,.. 



E. van Beveren, Z. 

arXiv:hep-ph/0602247 . 
P. KroU, Mod. Phys 

arXiv:hep-ph/0509031 ; 
oration]. Int. J. Mod 
|arXiv:hep-ex/0411081 ; 

F. Kleefeld, and G ~ 
B. Aubert 



* Electronic address: eefQteor. hs.uc.pt 
^ Electronic address: ,george@ist.utl.pt| 



Phys. C 21, 291 [1984) 

Lett. A 20, 2667 (2005) 
S. E. MiiUer [KLOE Collab- 
Phys. A 20, 1888 (2005) 
also see M. D. Scadron, 
Rupp, arXiv:hep-ph/0601196 

[BABAR Collaboration], 
"arXiv:hep-ex /0607082 

P. Colangelo, F. De Fazio and S. Nicotri, 
arXiv:hep-ph/0607245 

F. E. Close, C. E. Thomas, O. Lakhina and E. S. Swan- 
son, arXiv:hep-ph/0608139 

B. Zhang, X. Liu, W. Z. Deng and S. L. Zhu, 
arXiv:hep-ph/0609013 

In the one-PP-channel calculations of Ref. ^ and 
Ref. m (3rd paper), the D:o'(2850) pole was estimated 
at 2923 - i57 MeV and 2928 - i20 MeV, respectively. 
Here, with a restriction to the three PP channels only, 
the pole comes out at 2804 - i59.9 MeV. 



